Scattering and Recombination of Two Triplet Excitons in polymer light-emitting 

diodes 
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The scattering and recombination processes between two triplet excitons in conjugated polymers 
are investigated by using a nonadiabatic evolution method, based on an extended Su-Schrieffer- 
Heeger model including interchain interactions. Due to the interchain coupling, the electron and/or 
hole in the two triplet excitons can exchange. The results show that the recombination induces the 
formation of singlet excitons, excited polarons and biexcitons. Moreover, we also find the yields of 
these products, which can contribute to the emission, increase with the interchain coupling strength, 
in good agreement with results from experiments. 
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In recent years polymer semiconductors have made 
dramatic advances in light-emitting device performance, 
e.g., polymer light-emitting diodes (PLEDs) In these 
devices, following charge carrier recombination, singlet 
and triplet excitons are subsequently created. Compared 
to singlet excitons, which can decay radiatively and con- 
tribute to the emission, triplet excitons are nonlumines- 
cent. Although the precise ratio of singlets to triplets 
is currently debated, it has been established that triplet 
excitons constitute the dominant species in PLEDs (HI]- 
Therefore, research into the dynamics of triplet excitons 
in conjugated polymers is urgently required. 

Due to their long lifetime, a large number of the triplet 
excitons will accumulate in working PLEDs and have a 
significant probability of interacting and recombine with 
other particles. In our previous work pj, the recom- 
bining of a triplet exciton and a polaron was carefully 
investigated. Our results showed that the nonemissive 
triplet exciton is converted to emissive products, i.e., a 
singlet exciton and a excited polaron state. Therefore, 
the quantum efficiency of PLEDs can be enhanced by 
polaron-exciton recombination. In addition to this pro- 
cess, there is also another important case affecting the 
working characteristics of PLEDs, namely, triplet-triplet 
annihilation (TTA) fsl— fl~o|j - These works have demon- 
strated that in PLEDs delayed fluorescence (DF) orig- 
inates from the production of singlet excitons through 
the TTA process. Furthermore, Partee et al. 11] found 
that the DF intensity in films is higher than that in so- 
lutions. Accordingly, the triplet lifetimes are seen to de- 
crease from solutions to films. In a recently performed 
experiment, Ribierre et al. [12] studied the effect of in- 
termolecular spacing on the physics of exciton diffusion 
and light emission. They found the rate of TTA which 
is correlated to exciton diffusion, decreases with the in- 
termolecular distance. These two works llj, ll2| showed 
that both T-T annihilation by intcrmolecular recombi- 
nation and the total emission are higher in the presence 



of strong interchain interactions. Although some possi- 
ble mechanisms for TTA processes have been proposed, 
no detailed theoretical investigations of how triplet exci- 
ton recombination affects quantum efficiency have been 
reported. 

The goal of the present work is to investigate the re- 
combination of two triplet excitons on two coupled chains 
and to address the branching ratio of the products. We 
use a nonadiabatic evolution method at the unrestricted 
Hartree-Fock level within an extended Su-Schricffer- 
Heeger model including interchain interactions, with the 
extended Hubbard model for electron-electron interac- 
tions. Due to the interchain coupling, the electron and/or 
hole in the two triplet excitons can transfer between the 
two chains. As a result, the two triplet excitons can re- 
combine to form singlet excitons, as well as excited po- 
larons and biexciton states. These products are all lu- 
minescent due to radiative decay, which can enhance the 
efficiency of PLEDs. Furthermore, when compared with 
the results of experiments [8l-lll| , which indicate that the 
DF derives from singlet excitons formed in TTA process, 
our results show that the primary source of DF is excited 
polarons. 

The one chain model Hamiltonian we use is the Su- 
Schrieffer-Heeger (SSH) model 13[ with a Brazoskii- 
Kirova-type symmetry-breaking term [lij , including the 
effect of the intrachain electron-electron (e-e) interaction 
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Here, j = 1,2 is the chain index. The quantity t n is 
given by t„ = t - a(uj >n+ i - u hn ) + (—l) n t e with t be- 
ing the transfer integral of 7r-electrons in a regular lattice, 
a the electron-lattice coupling constant, and u^ n the lat- 
tice displacement of the atom at the n-th site along the 
j-th chain from its equidistant position. The quantity 
t e is introduced to lift the ground-state degeneracy for 
nondegenerate polymers. The operator c t n s (c J ,„ iS ) cre- 
ates (annihilates) a 7r-electron at the n-th site with spin 
s along the j-th chain; K is the elastic constant due to 
the a bonds, and M is the mass of a CH group. U gives 
the strength of the on-site Coulomb interactions and V 
gives the strength of the nearest-neighbor interactions. 

The two chains are coupled by the interchain interac- 
tions H' 0, 

H> = _ £*-L ( C l,n.s C 2,n, S +4,n, S C l,",s) 

(2) 

where t± stands for the transfer integral between sites 
labelled by the same index n on the two chains and V± 
is the e-e interaction term. In the following calculations 
V± is fixed at O.leV. 

The temporal evolution of the lattice is determined by 
the equation of motion for the atomic displacements, 

M U jtn = -K (2Uj >n (t) ~ Uj >n+ i(t) - Uj-, n _i(t)) 

+2a ]T [Pjnjn+l (*) " (*)] > ( 3 ) 

s 

where the density matrix p is defined as 

Pjn,j'n> (*) = E ®in,k (*) /fc,a*",»',ik(*). ( 4 ) 
fe 

with /fe jS being the time- independent distribution func- 
tion as determined by the initial electron occupation. 
The electronic wave functions $^ n k (i) are the solutions 
of the time-dependent Schrodinger equation 

(*) = E h fT.J' n '(*)*i',n',*(*)- ( 5 ) 

The coupled differential equations ([3]) and ([5]) can be 
solved with a Runge-Kutta method of order 8 with 
step-size control [17(. The parameters used here are, 
to = 2.5eV, a = 4.1eV/A, t e = 0.05eV, K = 21eV/A 2 , 
and M = 1349. 14eV • fs 2 /A 2 . 

In our numerical calculations, two 50-site parallel poly- 
mer chains opposite each other are considered: the CH- 
units are labeled as 1-50 on chain 1 and 51-100 on chain 
2. We start with the stationary solution to the system af- 
ter two triplet excitons are introduced, respectively, into 
chains 1 and 2. The initial bond configuration and elec- 
tronic structure can be obtained by minimizing the total 
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FIG. 1: (a) The staggered bond order parameters S n and 
schematic diagram of energy levels for the two triplet excitons; 
(b) The wave functions of the gap-state (from 49th energy 
level to 52th energy level). 

static energy of the two-chain system in the absence of 
the interchain interactions. 

Figure 1 (a) shows the initial staggered bond order 
parameters S n = (— l) n (tt n _i + it n +i — 2u„) /4 and a 
schematic diagram of the energy levels for the two triplet 
excitons. In polymers, the exciton is a composite particle 
consisting of an electron and a hole bound together by 
a lattice deformation, due to the strong electron-lattice 
interactions. With the two deformations, four localized 
energy levels emerge in the energy gap. We consider the 
two excitons to have opposite spins (with Sz = +1 and 
Sz = —1)- Obviously, there are two pairs of degenerate 
energy levels in the absence of interchain coupling. The 
corresponding wave functions are presented in Fig. 1 (b) 
[the 49th and 51st energy levels are localized on chain 1, 
and the 50th and 52nd energy levels are localized on the 
chain 2]. 

Then, turning on the interchain interactions, we sim- 
ulate the charge transfer process and the dynamic relax- 
ation of the recombination, with special attention given 
to the branching ratio of various species formed in the 
steady state. In order to avoid abrupt changes, the in- 
terchain coupling strength is smoothly increased over a 
period of 400fs, and then held at a constant value. As ex- 
pected, due to the interchain coupling, the electron and 
/or hole in the two excitons will exchange between the 
two chains. That is to say, a portion of the spin-up elec- 
tron and/or hole localized initially on chain 1 transfers 
to chain 2, and the spin-down electron and/or hole lo- 
calized initially on chain 2 transfers partially to chain 1. 
Because the total charges on the two chains are equal 
initially, the effects of the interchain interactions only re- 
distribute the electron and/or hole with different spin on 



3 



-X2 



singlet singlet 
(b) 

*r-B- — 

<^**y excited polaron polaron 

(e) 

*r-H- — zi 



triplet triplet 
(a) 



4^ 



negative positive 
bipolaron bipolaron 



(d) 



X 



-Zi 



biexciton ground 
(e) 

FIG. 2: The possible species of the recombination products 
of two triplet excitons. 



the two chains and do not change the quantity of charge 
on each chain. Correspondingly, the amplitudes of the 
two lattice distortions for excitons are unchanged in the 
dynamical evolution processes. When the interchain cou- 
pling strength eventually increases to its final value, the 
charge transfer also reaches a maximum. At this point, 
the system approaches a dynamically stable state. 

We now analyze the products of the recombination and 
determine the branching ratios for the formation of var- 
ious species in the steady state. The approximate eigen- 
states Xi (localized on one chain) and \i (localized on 
the other chain) can be separately obtained by ignoring 
the interchain coupling. It is found that such a structure 
is a mixture of sixteen possible states, i.e., an excited 
polaron localized on chain 1 and a polaron localized on 
chain 2, a polaron localized on chain 1 and an excited 
polaron localized on chain 2, etc. For simplicity, we di- 
vide these states into five species: (a) a triplet exciton 
on each chain, (b) a singlet exciton on each chain, (c) 
an excited polaron and a polaron, (d) a negative bipo- 
laron and a positive bipolaron, and (e) a biexciton and 
a ground state, see Fig. 2 [only the gap-states, Xi-> X\ 
for the particle localized on one chain, and X2 > X2 f° r the 
particle localized on another chain, are shown]. 

In order to ascertain the yields of these states, after 
each evolution step we can project the evolution wave 
functions |$ (t)) onto the space of the eigenstates of the 
system. The relative yield Ik if) for a given eigenstate 
\^k) is then obtained from 



I K {t) = |(*jc|$(t))| 



(6) 



The evolved wave function of the whole system |$ (t)) 



can be constructed using the single electron evolution 
wave function n k (t)) written as a Slater determinant. 

k) denotes the states of products shown in Fig. 2 (a)- 
(e). As an example, Fig. 3 shows how the yields of five 
species vary with time for the case t± = O.leV. It may be 
seen that the yield of triplet excitons decreases, and si- 
multaneously, the yields of the other four states increase 
as the interchain interactions are slowly increased. That 
is to say, the initial triplet excitons, due to interchain 
charge transfer, are transformed little by little into the 
other four states. After the interchain coupling strength 
reaches a constant value (at 400fs), the yields tend to be 
steady with only slight oscillations. In detail, the yields 
of the triplet exciton, the singlet exciton, the excited po- 
laron, the bipolaron, and the biexciton are about 66.8%, 
1.4%, 28.2%, 1.4%, and 1.4%, respectively. 

It should be stressed that the excited polaron is lu- 
minescent as is the singlet exciton because they have 
comparable transition dipole moments (about 30eA for 
the singlet exciton and 54eA for the excited polaron). 
Furthermore, the biexciton state, in which two electrons 
and two holes are trapped together by a lattice distor- 
tion, can decay to the ground state by emitting two pho- 
tons. Therefore, our results indicate that two nonemis- 
sive triplet excitons can be transformed into the emis- 
sive singlet exciton, the emissive excited polaron, and 
the emissive biexciton. Based on simple spin statistics 
arguments, the yield of the triplet exciton should be 75% 
and that of the singlet exciton should be 25% in PLEDs. 
When the recombination of two triplet excitons is taken 
into account, as stated above, for the case t± = O.leV, 
two triplet excitons have a probability of about 31% of 
being converted into emissive species (28.2% to the ex- 
cited polaron, 1.4% to the singlet exciton, and 1.4% to 
the biexciton, see fig. 3). Other processes also play a 
role. For example, in the recombination of the two ex- 
citons with parallel spin (Sz = +1 or Sz = — 1), new 
products do not form. In contrast, for the recombina- 
tion of the two excitons with Sz = +1 and Sz = 0, or 
Sz = — 1 and Sz = 0, or both Sz — 0, the emissive ex- 
cited polaron which is the primary product can be formed 
by charge exchange in the two chains. To summarize the 
above, in most processes where two triplet excitons re- 
combine (with Sz — +1,0, and —1), emissive products 
can be formed through interchain interactions. There- 
fore, the quantum efficiency of light-emitting diodes is 
enhanced, and the primary source of the increased ef- 
ficiency is the excited polaron. Consequently, the ex- 
perimental fact that the quantum efficiency exceeds 25% 
[18l [l9| can be understood on the basis of recombination 
of two excitons. 

With increasing interchain coupling strength, charge 
transfer between the two chains becomes easier, i.e., the 
transformation of various electron/hole states becomes 
easier. As a result, the products of the triplet-triplet re- 
combination, i.e., the excited polaron, the singlet exciton, 
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FIG. 3: The yields of five species vary with time in the case of t± = O.leV. 
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FIG. 4: The five species included in the final state as functions 
of the interchain coupling strength. 



the bipolaron, and the biexciton increase with increas- 
ing interchain coupling strength, while the yield of the 
triplet exciton decreases, see Fig. 4. This implies that the 
interchain interaction greatly increases the formation of 
the emissive excited polarons, singlet excitons and biex- 
citons, and thus the quantum efficiency of devices. Our 
calculated results are in goo d agreement with the exper- 



imental observations [111 Il2j. In ref 11 the delayed fluo- 



rescent of poly(p-phenylene vinylene) (PPV) and poly(p- 
phenylene ethynylene) (PPE) derivative solids and solu- 
tions was described. The results show that triplet-triplet 
bimolecular recombination induces the DF. Furthermore, 
the DF intensity in films is higher than that in solu- 
tions and the triplet lifetimes decrease from solutions to 
films. Reference 12 described the dynamics of triplet ex- 
citon quenching in phosphorescent dendrimers, in which 
the spacing of the emitting core was controllable. They 
found that as the distance between the cores increased, 
the rate of TTA decreased. These results showed that 
both T-T annihilation by intermolecular recombination 
and the total emission are higher in strongly interchain 
interactions. Therefore, our calculated results are indeed 
consistent with those found experimentally. 

Finally, the effects of the e-e interactions on recombi- 
nation are discussed. Depicted in Fig. 5 are the yields of 
the singlet exciton, which is equal to that of biexciton, 
and the excited polaron as functions of the on-site Hub- 
bard energy U for different nearest neighbor Coulomb 



interaction strengths V with t± = O.leV. The value of 
U is allowed to range from to 3.0 eV, and the relation- 
ships between U and V are U=3.0V (solid line), U=2.0V 
(dashed line), U=1.5V (dotted line). It is found that the 
emissive singlet exciton, biexciton and the excited po- 
laron are suppressed by U and V. The reason is that the 
binding energy of the electron and hole in a self-trapping 
excited state increases with the e-e interactions, making 
charge transfer between the two chains more difficult. 
Therefore, strong e-e interactions are unfavorable for the 
formation of the emissive species. 

In conclusion, we have simulated the two triplet exci- 
ton recombination processes using a nonadiabatic evolu- 
tion method. We have identified four types of products, 
of which the singlet exciton, the excited polaron, and the 
biexciton can decay radiatively. Our results indicate that 
the quantum efficiency of polymer light-emitting diodes 
is enhanced and can exceed 25% (the simple statistical 
limit) by the recombination of two triplet excitons. Fur- 
thermore, interchain interaction greatly increases the for- 
mation of emissive species, in good agreement with exper- 
imental observations. The effects of the e-e interactions 
on the recombination processes are also discussed. 
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